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The use of variational method in functional integral approach is discussed for fermion and 
boson systems with Coulomb interaction. The formal general expression of thermodynamic 
potential is obtained by Feynman path integral technique and representation of Coulomb in- 
teraction with functional integrals. Introduced additional complex field show to transform the 
problem to calculation of functional integrals containing third order vertices. The thermody- 
namic potential can be found from variational principle with respect to field cumulants. The 
calculation of the equation of state and critical properties is demonstrated for symmetrical 
plasma by variation of finite number of parameters in the propagator. 



I. INTRODUCTION 

The approach of functional integrals is used more 
and more frequently studying the properties of non- 
relativistic Coulomb plasmas. The Feynmann path inte- 
gral technique allows obtaining the properties of plasmas 
in systematic manner without utilising such approxima- 
tion as chemical picture or Pade formulaes. The path 
integral method takes the advantages of Green function 
formalism. Such ab initio method as restricted path inte- 
gral Monte Carlo simulations [Q provide plausible results 
at high temperature, where the experiments allows only 
some implicit measurements. The further development 
of these simulations is prospective for higher number of 
particles. The direct path integral method by Filinov ||] 
rigorously includes anti-symmetrisation but uses the ef- 
fective pair potential. Both methods are well suited for 
high temperature calculations. However, these simula- 
tions are still time consuming for complex fermionic sys- 
tems. Therefore, the deeper understanding of path inte- 
gral approach should be reached. The Hubbard-Schofield 
transformation has been a great leap forward for func- 
tional methods S. Recently, Brown and YafFe M have 
reported that additional integral over complex field is 
commendable calculating the action of a charged parti- 
cle. However, the approach in principe cannot be used for 
highly degenerate plasmas. The paper will show how the 
anti-symmetrisation effects could be included accurately 



by introduction of additional integration over complex 
field for large canonical ensemble. Tough, the mathe- 
matical methods to calculate the obtained formal general 
expression of partition function containing functional in- 
tegrals are insufficiently effective. There exists a varia- 
tional approach using the cumulative averages of fields for 
system in thermodynamic equilibrium. A universal func- 
tion of field cumulants that does not depend on physical 
model should be known for this purpose. The compar- 
atively inaccurate expression of the required function is 
obtained using the simplest diagram in expansion. The 
variational method is well suited for qualitative study 
of complex plasmas in the same manner as " Gaussian" 
packet in quantum mechanics choosing the appropriate 
class of variable functions. 



II. ANTI-SYMMETRISATION 

The Coulomb interaction belongs to a class with pair 
and positively defined potentials 



V(x) 



1 



e lkx , V k > 0. 



(1) 



For such kind of potentials, the Hubbard-Schofield trans- 
formation allows to write the density matrix of the sys- 
tem of n non-identical particles as || 



p n (X,X') 



2?</?(x, u) exp 



1 r u 



2V k 



exp 



uv(Q) ^ 2 
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Hpste^U), (2) 



(x\p j (U)\x')=p j (x,it',U) 



Dx(u) exp 



du 



mjX 2 (u) 



(3) 
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where integration is carried out over all configurations 
of the real field <p(x, u) in 4-dimensional space x € V, 
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u G [0, U = HP]; X(u) = {xi(u),x 2 (m), . . . ,x„(w)}; % - a 
charge of particle; tp^ is the Fourier component of electric 
field 



p(x) 



(4) 



The integration over real field y(x, u) does not depend on 
the number of particles, their charge and type of symme- 
try, Therefore, matrix pj(x, x', U) can be considered for- 
mally as density matrix of noninteracting particles placed 
in an external field iip(x.,u). The Coulomb interaction 
has been taken into account by integration over all con- 
figurations of this field. The result is similar to the Feyn- 
man interpretation of quantum electrodynamics by path 
integrals || , except the exclusion of the self- interaction 
part V(0) in (0). If some of the particles are identical, 
then either symmetry or anti-symmetry must be taken 
into account by summation over all permutations. The 
symmetrisation procedure should be performed for the 
system of noninteracting particles, according to the state- 
ment above. We are interested in the equilibrium prop- 
erties of the system of charged particles, such as equation 
of state or phase diagram. After symmetrisation Ij, the 
partition function of large canonical ensemble in path in- 
tegral representation becomes 



Z= fvip(x,u)exp ~T^2 I 

J i, J 



.£C/Trln(i-Ge^r 



Pj(U)) 



(5) 



where (x|l|x') = <5(x — x'); j now is an index of particle 
species (including different spin orientations); Q = ±1 
for bosons or fermions, respectively; /i*- ot - the chemical 
potential of j-type charged particle. The Fourier compo- 
nent Vfc of modified Coulomb potential V(x) is 



V(x) = -(1 



4?r 



hp 1 -j- ^2 



(6) 



where a large parameter T is introduced to avoid from 
the singularity of Coulomb potential at small distances. 
The constant self-interaction part can be added to 
chemical potential 



u tot 



(7) 



where \ij is the conventional chemical potential of the 
system. At the limit T — > oo, the terms proportional 
to T in the partition function must contract. The short 
distance divergences can be excluded also by introduction 
of fractional dimension [B . 



III. VARIATIONAL PRINCIPLE FOR 
CUMULANTS 

In the following section, the basic tool will be con- 
structed that helps to find the thermodynamic potential 
from variational principle for cumulative averages (cu- 
mulants) of some field. It is simply the method how to 
calculate the non-Gaussian type integrals. The idea was 
applied for the study of the phase transitions in Landau- 
Ginzburg theory Consider the partition function as 
a one-dimensional integral 



Z= d<pe Q{ip) . 
The statistical average of (^-function is 



dfS(f - ip)e<*W 



(8) 



(9) 



Therefore, the term Q(ip) which depends on the physical 
model can be isolated in partition function 



F = hxZ = Q(<p)-hi.(6{f-<p))f. 



(10) 



The partition function cannot depend on variable (p. 
Hence, it equals to its average: 



F = (Q(<p)) v -Qji(6(f-<p)) f ) 



(11) 



We can expand the average of ^-function in terms of cu- 
mulants using the integral representation of <5-function 



(*(/)>/ 



1 

2^ 



1 

2^ 



dh exp 



dh(exp[ihf]) f 



(12) 



where Kl = (f n ) c is the n-th order cumulant of the 
variable /. The relation to the usual average (/"), 



exp 



n=0 



(ihy 



,n=l 



Kl 



P, 



(13) 



is utilised in (|l2|) . Let us denote the n-th order cumulant 
of variable <p as K% — {<p n ) c - The variation of function 
F with respect to the total cumulant K n consists of two 
parts 



OF dF 



dF 



dK n dK[ dKl 



dF 



(14) 



= 0. 



It follows from the F independence of ip that dRV 
The derivative in respect to the cumulant K£, in compli- 
ance with (0, llj, fij) is zero, too: 



OF 
dKl 



dF 
dKl 
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2*(5(f-<p)) f 



— f°° dh ° P 



dh 



dP 



-ihcp 



(15) 



dP 

dh 7 5{h) = 0. 

dKl 



Subsequently, the variational principle is valid for ther- 
modynamic potential with respect to cumulants K n . The 
average (Q (</?)) ^ depends only on physical model and 
usually includes a finite number of K n . On the other 
hand, the average of the second term in ( |TT|) expressed 
by cumulants does not depend on interaction parame- 
ters, that makes the analogy of this term with an entropy. 
However, such universal expression is hardly obtainable. 
The expansion in terms of cumulants will be shown in 
next section using a diagram technique. 
If the partition function is the functional integral of some 
field tp(x) in a box x£F then it contains a product of 
^-functions 



{QM)h 




(16) 



in Fourier representation. Henceforth, the partition func- 
tion involves mixed cumulative averages (yki<^k 2 ■ ■ -)ip of 
physically different modes. But for a homogeneous sys- 
tem the part of interaction Q[tp] does not include mixed 
cumulants of ip^. Thus, one root after variation is always 
zero. 



IV. DIAGRAM EXPANSION ON FIELD 
CUMULANTS 

The diagram technique is discussed in a number of pa- 
pers both for classical plasma and quantum plasmas, e.g., 
H , [jl0| . The expansion is usually made in terms of De- 
bye screening radius. Such an approach is commendable 
for strict expansion at low densities. However, we need a 
specific diagram expansion on the cumulants in order to 
apply the variational principle. In first approximation, 
we will use the expansion in diagrams up to the square 
of charge. The final thermodynamic potential will like- 
wise contain the charge of particle up to its square. Of 
course, we will obtain similar results as in mentioned pa- 
pers though with different mathematical approach. Vari- 
ational principle provides essential advantage in compar- 
ison with virial or activity expansion studying the phase 
transitions in plasma. To apply the variational principle 
it is necessary to express the last term of ( |rT| ) in cumu- 
lants of some field tp. A general expansion exists probably 
for this quite universal mathematical task, but here the 
diagram technique is used for this purpose. The integral 
representation of the average of <5-function is (see (|l2|) 



(«(/-¥>)>/ 
1 

2^ 



= — / dh exp 



E 
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K n — ihip 



(17) 



where K n = (f n ) c = (</?™)c is the n-th order cumulant. 
Let us denote the cumulant K2 of the Gaussian part as 
G. The logarithm of @ can be formally expressed by 
means of diagrams as follows 



In <*(/-¥>)>,= In 



2tt 



dhexp 



h 2 G ^ (ih) n „ ., 



n^2 



E 




(18) 



where the first two terms follow from the Gaussian ap- 
proximation; the first graph represents the sum of all 
joint diagrams that does not contain ip, but the second 
graph all joint diagrams containing ip. We need to find 
an average of ( |l8| ) with respect to ip. The average (ip n ) 
can be written by cumulants in accordance with Vick's 
theorem 



(19) 



Therefore, the averaging of ( [l8| ) couples the first order 
vertices x to produce higher order vertices. The result- 
ing diagrams may cancel with the similar diagrams in the 



first graph of (|18|). Only certain kind of diagrams sur- 
vives. Thus, all diagrams that contain first order vertices 
or vertices joining vertex to itself (loop) cancels. As a re- 
sult, the average of (18) does not contain the cumulant 



K\ = (cp). The simplest remaining diagram consists of 
two K3 vertices: 



K 2 
2G 3 ' 



(20) 



All three edges cannot be crossed simultaneously, since 
the diagram would split apart before averaging with re- 
spect to ip. Therefore, the first few terms in expansion 



4 



(21) 



The intrinsic energy part in Fourier representation usu- 
ally does not contain cumulative averages of order higher 
than second. Then, the extremum condition 



d In z 

dK 3 





yields A3 — 0. However, that could probably violate for 
additional extremums if the next diagrams are consid- 
ered. 



V. COMPLEX FIELD FOR FERMIONS 

The symmetry properties of particles are accounted 
in the partition function via the logarithmic term in 
(||). It seems that such a construction spoils its solv- 
ability due to the slow convergence of Taylor series for 
strongly coupled plasma. Nevertheless, let us remember 
that any nonsingular quadratic matrix satisfies equation 
TrlnA = lndet A. At the same time, the determinant 
of matrix A can be represented as the path integral over 
all configurations of complex field ^(x) [O: 



1 



det A 



(22) 



2?W(x) exp 



- / dx 



J dx'^(x)K(x,x')^*(x') 



This leads to the idea that the logarithmic part of the 
partition function (|^) for fermions can be represented as 
a functional integral 



exp[— Tr ln(l 



pAu))\ 



= jv^xjexp -Jdx\^(x)\ 2 (23) 
+e P4 at +in [ dx J dx'* j {x)p j (x,x' ) U)**{x') 

The multiplier e lTr is added in ( p3| ) to change the sign 
for fermions. The change of the sign is unnecessary for 
bosons: the difference will be shown latter. Note an anal- 
ogy of function VPj (x) with the stationary wave function 
in quantum mechanics. As mentioned above, the inte- 
gration over complex variables \& is introduced also in 
H for canonical ensemble. After splitting the imaginary- 
time interval in JV > 1 small parts e = the kernel 
/9j(x, x',U) (3) of linear integral operator pj becomes a 
convolution [6] 



Pj(x,x',U) = J dxL N -i . . . J dx 2 J dxipj(x,x N -i,e)...pj(x 2 ,Xi,E)pj(xi,x',e). 

Henceforth, the logarithmic part ( p3] ) of the partition function can be transformed to a Gaussian type functional 
integral over all configurations of complex field *ifj (x, u) in 4-dimensional space Jll| 



J dx J dx' i Sfj(x,u-\-e)exp 



U 



P] (x,x',e)^*(x',u) , 
u = 0, e,..., (N -l)e. 



For Coulomb system with density matrix pj(x, x', e) (g), 
the charged particle can be considered as placed in imag- 
inary field iip, their Hamiltonian being Hj — —^-A x + 
iqjtp (a non-Hermitian operator). The statistical density 

matrix P j(x,x',u) satisfies the equation = —HjPj 
analogous to Schrodinger equation. The expansion of this 
density matrix and small e is P j(x, x',e) for small e is 



-— q y( x >,u))S(x-x>). (24) 

Hence, the partition function (|^) becomes a functional 
integral over all configurations of the real field <^(x, u) 
and the complex field \& j (x, u) 



Pj(x,x',s) w 1 



eh 
2m i 



: %^(x',«) 
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Z = / XV(x, ■ 



exp 


e |y k (ti)| 2 
ft Z^k,u 2Vfc 




n 3 -/2?* i (x,u)exp 




)-f*,(x,U + £)W J -]^(x, u) 



(25) 



2m, 



I 



ink 
IT' 



(26) 



The number of complex fields coincides with the number 
of particle species. The Fourier transformation of these 
fields is useful: 



9 
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IT' 



n = -M,-M + 1,...,M 



(27) 



assuming for convenience that JV = 2M+1 is an odd inte- 
ger. The e term in operator Hj ( p6| ) should be neglected, 



J 



but for systems with Coulomb interaction it survives and 
accounts for self-interaction. As we shall see in (|5l]), the 
average Gk, w = (|<^k,w| 2 ) at |w| — > oo equal to the inter- 
action potential Vfe. Thereafter, the self-interaction part 
arises, according to the Fourier expansion of Coulomb 
interaction (^J): 



fN 

lime<^(x,u)) = lim — Y,V k = hV(0). 



(28) 



The self-interaction part cancels if the chemical potential 
of Coulomb system u* ot , (0), is subtracted in (p6|). The 
partition function (|2q) for small e becomes 



Z = I T><p(x, u)- 



exp 


\<Pk,u>\ 2 

Z^k,w 214 




n i /X>*i(x,u)exp 


- E k , w (^k,.|^ )k) .| 2 + v ix * iiq)X *5, k+qiW+x )" 



(29) 



where 



oy,k,< 



A^(e" 

id 



: /3 



1) — in - 

h 2 k 2 



(30) 



2m i 



The term — i7r is not present for the system of bosons. 
Moreover, the integration over complex field ^j(x,u) in 
partition function Z for bosons is present in numerator. 



VI. VARIATIONAL PRINCIPLE FOR 
FERMIONS 

According to the previous section, the following inte- 
gral is to be found in quantum statistics for Fermi sys- 
tems: 



dtp 



exp[Q(tp)] 



J d*exp[$(v5, |f| 2 )]' 



(31) 



where = dRe^! dlm^f ; ip is real variable, but '5 - com- 
plex. The average values of (5-function are 

W-<p)8(*-Q)) f , 9 



(«(/-¥>)>/ 



1 exp[Q(y)]exp[<l>(^|e| 2 )] 
^(/d^iexp^,!*!^)]) 2 

1 exp[Q(tp)} 



Z J dtfiexppfoltfil 2 )]" 



(32) 



(33) 



Therefore, the partition function can be split using 5- 
functions (|3^, |33| ) in analogy with simple field of variable 
p in (FT 



InZ 



(Q(^)) v -($(^|e| 2 )) ¥ ,e 
- 2(ln(5(f-tp)) f ) v 

+ (M8(f-<P)S(*-Q))t,9) v ,e 



The conclusions are similar to those made in sect. Ill, 
i.e., the variation of partition function with respect to 
cumulants (f l ^ m ^* n ) c must be zero. In our case, 
we should consider the fields in coordinate-temperature 
space, <£>(x, u), ^(x, u) as in (|l6|). The first two terms 
follow from functional integral representation (EP 



(Q[<p]) v -<<%, |6| 2 ])^ e = 



\\<Pk,t>>\ )<P 



2V k 



6 



+ E «^(|e,, k ^| 2 )e (34) 

One recognises here the terms representing the energy of 
electric field, the energy of noninteracting fermions, and 
the particle-field interaction characterised by the third 
order vertex (</?(x, u)\l/ J (x, u)5'*(x, u)) in the coordinate- 
temperature representation. 



VII. DIAGRAM EXPANSION ON CUMULANTS 
FOR FERMIONS 

The expansion of 

(\n(S(f - <p)8(9 - Q)6(** - 0*)) /,*,*.) ¥ , ) e,e«, (35) 

f,<?eR, *,6eC 

in terms of cumulants is necessary for fermion plasma. 
By using the same scheme as in the previous section, the 
integral representation of <5-function gives an expansion 
in cumulative averages 



1 



(27T) 
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dhdadb exp 



E 

ni+ri2+"3>l 



(i/i) ni (L4)™ 2 (M*)" 3 
ni!n 2 !n 3 ! 



(/m^n^*»3) c _ ihtp - iAS - iA*Q* 



where a = ReA and b = ImA are real integration vari- 
ables. Since W and \P* are complex- valued, the cumulants 
(J n nj/ n 2\]/*™3^ c are zerQ ^ n ^ ^jjj ers f rom n3 ^ Moreover, 

the average (ip) is zero for neutral plasma. Therefore, the 
first nonzero summands of the sum in previous equation 
are 



E 



n 1 +n 2 +n 3 >l 



(ih) ni (iA) n2 {iA*) r 

riling!?^! 



h 2 G 



- AA*n-ihAA*K + (36) 



where II = (\I r \f r *) c is the propagator. K = (/'& 1 E r *) c is 
the vertex representing the particle-field interaction, 




where the wiggled edge corresponds to electric field /, 
and edges with arrows to complex fields W and . 
Higher order cumulants does not survive after variation 
as the interaction part in the partition function (|3j) does 
not include them. The presence of K vertex can be in- 
cluded by means of diagrams. All diagrams that contain 
either first order vertices or loops cancels. The simplest 
remaining diagram contains two K vertices: 



50 



K 2 



2Gnn 



(37) 



Ignoring higher order diagrams, the average becomes 



(\n(6(f - <p)5(* - Q)5(V* - &*)) fw) v , e ,e* 



= --(ln(27rG) + l) 



(ln(-i27rn) + 1) 



K 2 



2Gnn 



(38) 



The next diagrams containing four vertices K are not 
included. The accurate integration of the single integral 
(|3lj), containing only mentioned cumulants G, II, and 
K in the interaction part, shows that the sum of all di- 
agrams gives I hi f 1 + fjnn) > which agrees with ( |38| ) at 
small K . However, the interpretation of the result is non- 
trivial for many-dimensional integral without the use of 
diagram technique. K is approximately proportional to 
the charge of a particle qj. Therefore, the higher or- 
der diagrams as fl37| ) can be neglected considering the 
Coulomb interaction as a small perturbation. Therefore, 
the partition function of fermion plasma is 



InZ = const H — ( m 



2f- V V k 

K.,UJ 



E ( mn j,k.UJ + 1 - ";.k... H,.k.. 

w 



(39) 



E 



K 2 



j,k,w,q,x 



2Gk,ajnj 5 q,^nj" ; l c +q,aj+x 



where Gk, w = (|</5k,w| 2 ) is cumulant of electric field, 



propagator, and i£j - ,k,w,q,x = 
cumulant of field-par- 



n i)k , w = (|* Jlk ,u,r) - 

(^k,u,*j,q, X! *j,k+q,u,+ X 

ticle interaction. The constant can be readily deter- 
mined by comparison with an ideal system, where In Z = 
J] j k ln[l + e _c 3 k ]. The first sum in the partition func- 
tion corresponds to contribution of field, the second one 
to noninteracting fermions, while the last ones are due 
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A 



FIG. 1: Summation over u> 

to fermion- field interaction. The variation of -ftT^k.u.q.x 
gives 



K 



V 



Gk.o,n,- 



<rij,k+q.w+x 



(40) 



Thus, the diagram in ( |37| ) is the only one which is pro- 
portional to the square of the charge assuring that self- 
interaction part contracts exactly in sect. The func- 
tions Gk, w and rij k,^ should be obtained from the mini- 
mum of the thermodynamic potential, too. 



VIII. APPROXIMATION OF FERMION 
PROPAGATOR 

The correlation function Gk, w is real but propagator 
IIj.k,w - complex-valued. It seems that the main difficul- 
ties are connected with nj.k,w Therefore, we will ap- 
proximate the propagator so that it depends on to in the 
same way as reciprocal of a^k,^ does: 



n 



1 



N{e lule — 1) — ztt + • 



(41) 



where the real function Cjjk does not depend on lu. If wc 
try to give some physical meaning for this function, it can 
be considered as self-energy. The variational principle is 
valid also for the new function Cj.k, which may be further 
simplified. The density of j-type particles is 



1 <91nZ 
V/3 dfij 



1 v 1 ^ 



"i,k, (42) 



The summation can be represented by a set of N (see 
(p7|) ) points distributed along a circle as shown in Fig. [|. 
Let us divide the circle in two arcs: A, where \we\ <C 1; 
and B. Discrete summation is necessary in the arc A. 
The arc includes infinite number of points at N — + oo. If 
the angle of this arc is sufficiently small, the propagator 
is approximately Ilj^k.w ~ (*7r(2?i~l)+c : j,k) _1 . In the arc 
B, the summation can be safely replaced by integration 
over cj) = we. The integration limits are — tt and 7r, when 
the angle of the arc A is small. The propagator in the 
arc B is ILj t k,v ~ (-/V(e ,( ^ — Hence, one obtains the 

well known form of Fermi-Dirac distribution: 



n,-k = -r- 



1 

2n 
1 



i<p 



oo 

E 

71= — OO 



1 



Z7r(2n — 1) + Cj.k 



1 + e c i.' 



It is a consequence of approximation ( pl| ) 
tion of (E3) yields the sum of logarithm: 



(44) 

The integra- 



E 



inn 



i,k,w 



ln[l 



const. 



(45) 



The last sum of the partition function ( p9|) contains the 
product of the propagators. The single sum of the prod- 
uct at small M is 



n^q^II^q+k, 



n 



j.q+k 



i2im + c 



i.q+k 



(46) 



\n\ < N, 



because the integral along the arc B (see Fig. |l|) cancels. 
The double sum is a consequence of (43): 



n ? ,q,^rij,q+k,x+^ — n j,q n j,q+k- 



(47) 



Thereafter, the partition function essentially simplifies 



'i,k - E n ^k.- - E N f e iue _ 1) _ i7r + c . k • 

(43) 
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G*k,n 

v k 



J2 (Ml + ^ k ] + (c,, k ~ c^Kk) 
J,k 



E 9 i ^ ( " <i"/ <i-k - E 

j,q,k \ n 
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0k,; 



j,q+k 



i27rn + Cj,k+ q 



(48) 
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where the summation index n now runs from — oo to oo. 
The last term does not contain singularities at n = 0, 
because nj^, (|44|), monotonously depends on c^k. The 
partition function does not contain divergent terms and 
so it is not necessary to introduce screened potentials. In 
the second sum, one recognises the logarithm of thermo- 
dynamic probability 



ln[l 



= -n^khin^k - (1 -nj,k)ln(l -nj,k)- (49) 

The variation of the correlation function of electric field 
yields 



G 



k,n 



Vu 



G 
V 



k . n 



(50) 



2 ( n i:q+k n j,qJ ( c i,k+q 



V 2^ (2 7 rn)2 + ( Cj . k+q -c i , q ) 2 



3,1) 



This relation proves the limit Umi^i—Kx, Gk. u — Vk nec- 
essary for exclusion of self- interaction in (|2S|). Subse- 
quently, the partition function depends only on one func- 
tion Cj.k: 

InZ = ~ (5 k|Tt - ln[l + 5 k , n ]) 

k,n 

+ (M 1 + e~ Ci ' u ] + (Cj,k - 4f k )^,k) (51) 

The variation of c-^k seems more complicated. Hence, 
the restriction of c^k in a certain class of real functions is 
necessary. This function may not be smooth at zero tem- 
perature, similarly, as the spectrum of some substance 
becomes sharper decreasing the temperature. Note, that 
one can in principle obtain the ground state energies for 
a given molecule performing the limit T — > in the ther- 
modynamic potential. 

At low density limit, one obtains a classical Debye-Huckcl 
screening in position space: Gq(x) = e~ ax jx. The last 
term of ( |54| ) contains the product of distribution func- 
tions n^k and it is negligible at low densities. Hence, 
the minimum of thermodynamic potential at low density 
limit gives the difference 
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-KD, 



(52) 



where up is Dcbye-Hiickel screening radius. This differ- 
ence does not depend on the wave vector k. It suggests 
that the presence of Coulomb interaction in the system 
is accountable by scaling of the chemical potential. 



0.6-1 
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FIG. 2: Correlation function of electric field at T = 0.01 and 
various densities for symmetrical plasma. Atomic units. 



A. Boson propagator 

The function <Xj )kin , (|30|), for bosons differs by — in: 

N{e %2 7r - 1) + cf k . (53) 



oy,k, 



The approximation of type ( [ll| ) leads to the well known 
form of Bose-Einstein distribution: 



%,k 



2tt 



J_n e iz - 1 + i2nn + 



c 7 ,k e c J- k — 1 



The partition for bosons can be obtained similarly as 
for fermions keeping in mind the comments for boson 
particles in section [v|. Finally, one gets the partition 
function 



ln2 = 



ln[l + a k ,„]) 



c }fkK',k) 



3,k 

JL 

2V 



X! s;'/; l i"; .|";.<i k- 



Ci' 



(54) 



(55) 



where Q = ±1 is for bosons and fermions, respectively. 
This partition function can be applied for plasmas con- 
sisting of both fermions and bosons, e.g., deuterium 
plasma. 



IX. SYMMETRICAL PLASMA 



Usually, only the simplest Coulomb systems are of 
theoretical interest, e.g., one-component, electron-hole 
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FIG. 3: Pressure isotherms vs. chemical potential for sym- 
metrical plasma. The intersection point of curve corresponds 
to phase transition point. 



(symmetrical plasma), hydrogen, deuterium and helium 
plasmas. Here only the quantum symmetrical plasma 
is considered as an example, i.e., all species of particles 
are fermions. The symmetrical plasma is of interest be- 
cause the controversial discussion exists about the loca- 
tion of the critical point for its first order phase transi- 
tion both in classical model of hard cores ial and quan- 
tum case . The unpolarised symmetrical plasma has 
four species j of particles because of two possible spin 
orientations. Despite the propagator has already been 
approximated in (41), it is still difficult to obtain the 
function c,,k from the variation of thermodynamic po- 
tential (pTJ) . Therefore, only two parameters are variated 
for each specie of particle: 



(56) 



Thus, a scaled chemical potential ftj and an effective mass 
rhj = rrijrjj /rjj are chosen as variational parameters. The 
correlation function of electric field is calculated on the 
basis of ( pl| ) for every n. For convenience, the atomic 
units are used, and temperature is measured in the units 
of energy (hartrees). The symmetrical plasma can be ei- 
ther electron-positron plasma or more familiar electron- 
hole plasma. However, the number of particles in both 
systems does not conserve due to annihilation of antipar- 
ticles and recombination of electrons and holes in semi- 
conductor. The phase transition of electron-hole plasma 
is established both theoretically, e.g., in Il3], JlH] 

and experimentally [ pH) . 

It would be helpful to know what kind of the correlation 
function Gk, w , (pl|), one obtains for fermions, when Cj k is 
chosen according to (|56"|). The correlation function Gk,o 
is of interest corresponding to an average with respect to 
imaginary-time since only terms with lu — are present 
in classical plasma. Fig. || shows the correlation func- 
tion in position space Gq(x) at temperature T — 0.01 




0.01 
temperature 

FIG. 4: Phase diagram of symmetrical plasma. Variational 
parameters: (1) - jlj; (2) - flj,rhj. Atomic units. 



and different densities. The correlation function trans- 
fers from screened behaviour, Gq(x) ~ V(x) — e~ ax /x, 
to oscillating one, when the Wigner-Seitz radius is com- 
parable with the Bohr radius. The method does not yield 
particle-particle distribution functions. 
Since the partition function is linear with respect to the 
chemical potential /ij, it is better to fix the variational 
parameter ft while fij follows from the system of varia- 
tional equations. The effective mass monotonously in- 
creases increasing the density. Thus, the presence of 
weak Coulomb interaction smoothes the step-like Fermi- 
Dirac distribution function. The phase transition point 
can be found plotting the pressure vs. chemical poten- 
tial (see Fig. ||) at temperatures below critical one. The 
intersection-point of isotherm corresponds to the coex- 
istence of two phases in accordance with phase equilib- 
rium condition. Loop with negative pressure corresponds 
to unstable densities. The obtained phase diagram of 
symmetrical plasma is shown in Fig. |J. The critical 
temperature is considerably lower, when purely quantum 
terms Gk,n with n^O are neglected. Obtained critical 
parameters of non-annihilating electron-positron plasma 
are T c = 11500X, P c = 1.6 • 10 8 Pa, n c = 5.7 • 10 : 



27, 



(for curve (2)). The critical temperature appears to be 
slightly higher than the corresponding quantum statisti- 
cal result in T c — 7186X, while the density is ap- 
proximately the same. The hydrogen-like bound states 
are introduced directly in [fl3| , while the ground state en- 
ergy in the presented functional integral model appears 
implicitly with slightly different value at T — > 0, where 
the method does not work. One notices in Fig. [|that the 
equilibrium concentration of the dense phase (curve (2)) 
is inversely proportional to temperature. That would 
lead to collapse of symmetrical fermion plasma at zero 
temperature. However, the partition function ( pi"} ) is, of 
course, not valid for strongly coupled plasmas. Simula- 
tions made by restricted path integral method jlB) show 
that such phenomena as Bose condensation of excitons 
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and biexcitons takes place in dense electron-hole plasma. 



X. CONCLUSIONS 

The paper is devoted mainly for elaboration of varia- 
tional methodology that helps to investigate the quantum 
Coulomb systems with Feynman path integral technique. 
The principal points the proposed scheme are: 

1. The representation of the partition function of 
fermion gas by functional integrals is obtained us- 
ing Feynman path integral and formally introduced 
integration over complex field. This integral con- 
tains only third order vertices that simplifies fur- 
ther use of diagram technique. 

2. The form of the functional integral differs essen- 
tially for Fermi and Bose systems. For example, the 
integration over complex field for fermions stands in 
denominator, while for bosons - in numerator. Such 
a mathematically formal representation solves the 
problem of symmetrisation. 

3. The self-interaction part is excluded using modified 
Coulomb interaction potential at small distances 
and scaled chemical potential. The limit to accu- 
rate Coulomb interaction potential is performed in 
final expressions. 

4. The algorithm for expression of thermodynamic po- 
tential by average values cumulants of real and 
complex fields is elaborated. Those values of cu- 
mulants can be further found from the minimum of 
thermodynamic potential. 

5. The phase diagram and critical point are obtained 
using the simplest approximation of thermody- 
namic potential. For this reason, variation only 
of chemical potential and effective mass for corre- 
lation functions of particles is performed, while the 
minimisation in respect to cumulant of electric field 
is maintained exactly. 
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